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ABSTRACT 


In this paper, scalar limitations of diffractive optic components are investigated using coupled wave analyses. Results are 
presented for linear phase gratings and fanout devices. In addition, a parametric curve is given which correlates feature size 
with scalar performance. 


1. INTRODUCTION 

Diffractive Optical Elements (DOE) are utilized in optical computing architectures, fiberoptic systems, imaging cameras and 
numerous optoelectronic devices^. The major advantages are realized in coherent systems, due to the inherent chromatic 
behavior of these elements. Since DOE’s are fabricated using standard photolithographic equipment, the feature sizes can be 
on the order of wavelengths. However, when the features are a few wavelengths, as is the case for fast microlens arrays, the 
scalar assumptions are no longer valid. To address this issue, exact electromagnetic modeling must be incorporated in the 
design process. 

In this paper, the rigorous coupled wave formulation is implemented to determine where the separation is between the scalar 
theory and exact formulations. First, an analysis is presented for a discrete linear phase grating to demonstrate the 
polarization and coupling effects as a function of period size. Second, a Dammann fanout device is analyzed to demonstrate 
the effects on beam uniformity as a function of both polarization and period size. Third, an attempt is made to characterize 
the effect of feature sizes characteristic in multilevel structures. 

2. THEORY 

Diffractive optic components rely on phase perturbations in the transmitted wavefront to diffract the light Quite often this is 
performed using 2ir structures etched into a glass substrate using standard photolithographic techniques. Figure 1 illustrates a 
collapsed diffractive lens. This plate can be viewed as a series of zones diffracting energy into a prescribed diffraction order, as 
the zones move away from the center they become smaller since they are required to diffract energy through larger angles. 
Therefore, the smaller zones can be modeled as a dielectric grating since adjacent zone dimensions do not significantly change. 
For a multilevel diffractive optic, the -1 diffraction efficiency can be approximated using scalar diffraction theory as follows: 



( 1 ) 


where, T represents the transmission coefficient squared and N is the number of phase levels used to approximate the profile. 
To determine the diffraction angle for a given zone, the grating equation is employed 
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where X is the wavelength and n is the index of refraction. The period is the zone size which is modeled as a grating unit 
cell. Another side effect of the small zones is that the scalar, or paraxial, assumptions associated with diffraction efficiencies 
are no longer valid and more robust techniques are required. Therefore, Maxwell's equations must be implemented to 
accurately model the polarization sensitivities and resonances present in small structures. Although numerous formulations 
exist under diverse titles, the majority require solutions to a system of differential or integral equations. The most common 
technique used for dielectric grating profiles is the coupled-mode approach, wherein a system of differential equations is 
solved t 2 ]. The first step is to expand the Fields in the incident and transmitted regions in terms of a common basis set 
representational of the floquet modes 
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Since the modulated region is a multilayer structure of air/dielectric rectangles, it is expanded in the same basis set This 
yields the following set of equations for the Fields and modulated index: 
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operating on these field components with the curl equations, a set of coupled differential equations is formed for the modulated 
region ; — — . , , • - • • • • • . ■ - 
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Then, an eigenvalue/vector computation is performed to derive the solution set for each sublayer. 
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After this is accomplished, the fields are matched to the incident and transmitted regions for computing the Rn and Tn 
coefficients. These are then used to compute the diffraction efficiencies and field values at any desired point. 


3. NUMERICAL RESULTS 


3.1 LINEAR PHASE GRATING 

In order to determine where the scalar theory breaks down with a linear phase grating, a multistep grating : was “ Jg 
the coupled wave theory. Figure 2 illustrates an eight phase level structure with the polanzauonsdefined ^ E and H where 
E is parallel to the grooves and H is when the E lies in the plane of incidence. The computed diffraction efficiency for the -1 
order is illustrated in Figure 3 as a function of period size and index of refraction. As can be expected for periods 
approximately a few wavelengths in size, scalar theory breaks down; whereas, for periods greater than 20 wavelengths the 
performance is close to scalar. Moreover, as the index is increased the scalar approximation is realized at smaller periods. 

As was previously mentioned, the primary mechanism in the diffraction process is the phase modulation; therefore, 
to diffract the energy into the -1 order, the phase should be a stepwise Unear approximation across the unit cell To verify 
that the phase was Unear after immediately passing through the profile, a scan of the phase was computed for both large and 
smaU ooriod sizes Figures 4 and 5 represent the results for both polarization states at penod sizes of 2 and 20 wavelengths 
respectively. The smanJSeriod demonstrated a complex phase profile for both E and H modes which signified a departure from 
scalar theory. In the 20 wavelength case, both polarizations demonstrated a well defined linear phase function close to the 
scalar approximation. 


3.2 FANOUT DEVICES 


Another area of interest is in fanout phase grating structures. Although numerous techniques exist for beam splitting, one of 
the most common is a Dammann grating as iUustrated in Figure 6. Due to the thickness of the grating structure, the 
transmitted Ught assumes a 0 or JC phase shift which is then used to diffract the Ught into a predetermined number of diffracted 
orders with uniform intensity^. However, as the number of desired spots is increased, so is the number of transition points 
in the grating profile. The increase in the number of transition points forces smaller features m the grating design which 
must be investigated using rigorous modeling techniques. 


As an example, a five beam Dammann device was analyzed as a function of feature size. Figure 7 illustrates the unit cell 
which was designed in glass (n=1.5). This design was optimized using a steepest descent algorithm with scalar assumptions 
for a spot uniformity of less than 1% and an average DE of 14.9%. Uniformity is defined as follows: 


%Unif = 100x| Pm p~ Pmin -j 


(7) 


Figures 8 and 9 illustrate the spot uniformity and average DE as a function of period size for the 5 beam array. These results 
demonstrated polarization sensitivities and an unacceptable beam uniformity for period sizes up to 50 wavelengths. Therefore 
to obtain acceptable performance out of a Dammann device, the period must be hundreds of wavelengths when there are small 

features. 

3.3 FEATURE SIZE STUDY 

Many concerns exist as to when scalar theory actually breaks down and how do you design around this phenomena. 
Additional studies have implicated that it is not the period size, but the actual feature size relative to a wavelength that causes 
the separation from scalar. In an attempt to verify this effect, numerous analyses were performed on a biphase-level structure 
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with alternate fill factors as illustrated in Figure 10. Seven fill factors were analyzed: 0.9375, 0.875, 0.75, 0.5, 0.25, 0.125, 
and 0.0625. Since most of the scattering occurs at the comers of the structure, only the E mode was analyzed with the 
coupled wave theory. The -1 order diffraction efficiency was compared to Fourier theory as a ratio. Figure 11 illustrates the 
result for an index of refrection of 1.5. The implications of these results signify extreme fluctuations in the DE as a function 
of feature size and as the index of refraction is increased, the effect is diminished. 

4. CONCLUSIONS 

In this study, the rigorous coupled wave approach was used to investigate the scalar limitations of diffractive optics. This 
analysis verified that coupling effects influence the diffraction efficiencies of the multistep phase gratings. Additionally, 
polarization sensitivities were observed for Dammann gratings with small features! “Considering these findings, rigorous 
analyses must be performed when designing DOE’s with feature sizes on the order of a wavelength. Areas of additional 
studies include the dispersive effects of multliphase structures, and angle of incidence effects. Further developments in the 
analytical tools are also required to adequately predicts the three-dimensional diffraction of microlens structures. Moreover, 
alternate formulations from the coupled mode approach are needed due to the increased ode coupling of a three-dimensional 
structure. 
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Figure 1. Diffractive col lapsed-phase 
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Figure 6. Optical implementation of a Dammann grating. I 
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Figure 7. Dammann 5 beam unit cell. 
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Figure 8. Computed beam uniformity for a 5 beam Dammann device. 
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Figure 9. Computed average fo t a 5 beam Dammann device. 
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Figure 11. Exact and scalar comparisons for E mode as a function 

of feature size. 

(fill factors > 0.5 are given in the lower half) 




